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We carry out a theoretical study on the isotropic-nematic phase transition and phase separation in 
amyloid fibril solutions. Borrowing the thermodynamic model employed in the study of cylindrical 
micelles, we investigate the variations in the fibril length distribution and phase behavior with 
respect to changes in the protein concentration, fibril's rigidity, and binding energy. We then relate 
our theoretical findings to the nematic ordering experimentally observed in Hen Lysozyme fibril 
solution. 



PACS numbers: 87.14.em, 87.15.Cc, 05.20.Gg 



I. INTRODUCTION 

Amyloids are insoluble fibrous protein aggregations 
stabilized by a network of hydrogen bonds and hydropho- 
bic interactions 0,0,0,0). They are intimately related to 
many neurodegenerative diseases such as the Alzheimer's 
Disease, the Parkinson Disease and other prion diseases 
0. Furthermore, it has recently emerged that non- 
pathogenic amyloid fibrils possess great technological po- 
tential. In particular, amyloid fibrils have been employed 
as nanowire templates [6|, l7| , were shown to possess great 
tensile strength [8|,[9| and complex phase behavior similar 
to liquid crystals [loL ITU [l^ | . Given these extraordinary 
properties, it is highly desirable to investigate how one 
may exploit amyloid fibrils as functional materials. Here, 
we study theoretically the isotropic-nematic phase transi- 
tion in amyloid fibril solutions by combining the physics 
of self-assembled linear structures, as studied in cylin- 
drical micelles (see [H, [TiJ and the references therein) , 
and the physics of the nematic ordering in charged rods 
[T5L [l6| . We then apply the formalism to a specific ex- 
ample - Hen Lysozyme (HL) fibril solution, and discuss 
agreements between theory and the experimental results 
in 0. 

In the next section, we introduce a toy model for amy- 
loid fibrilization and review briefly the physics of nematic 
ordering in self-assembled rods. In Section III, we apply 
the theoretical formalism to Hen Lysozyme (HL) amy- 
loid fibrils and estimate all of the model parameters from 
previous experimental studies. We then discuss the limi- 
tations of and predictions from the model in Section IV. 



II. THE MODEL 

We assume that the monomers self-assemble into the 
fibrillar form through two different interactions: i) A- 
type interactions of strength a which are directed longi- 
tudinally along the fibrillar axis, and ii) B-type interac- 
tions of strength f3 which are lateral to the fibrillar axis. 



For amyloid fibrils, the A-type interactions would corre- 
spond to the hydrogen bonds among the beta strands and 
hydrophobic interactions between the side chains packed 
between the beta sheets; and the B-type interactions 
would correspond to the inter-cross-beta-sheet interac- 
tions (cf. Fig. H). We further assume that all fibrils 
are formed with the same number of filaments, 7, which 
is peptide specific and is constrained by the chiral nature 
of the cross-beta sheets [13, [3 • 

Given N monomers in a volume V of solution, we 
denote the fibrillar aggregate consisting of s monomers 
by N s , i.e., J2 S S ^ S = N. We consider only the fib- 
rillar species and ignore the free energy contributions 
from monomers and oligomers not in the fibrillar form. 
This assumption is satisfied if the concentration is much 
higher than the Critical Fibrillar Concentration (CFC) 
(so that the monomer concentration is negligible) , and if 
the CFC is lower than the critical concentrations of other 
oligomeric species [19( • These conditions are met if fibrils 
constitute the most dominant species in the solution. 

As each monomer is a peptide, there are intrinsic inter- 
nal degrees of freedom which contribute to the partition 
function. This is equivalent to the contribution of con- 



FIG. 1: (Color online) Schematic diagrams of a monomer 
and a fibril, (a) The monomers interact with each other via 
two types of directional interactions indicated by the blue 
(dark grey) arrows (j4-type interactions) and the red (shaded 
areas between the beads) patches (B-type interactions) (see 
text), (b) A fibril is formed by joining the blue arrows and red 
patches. The directionality of the A-type interactions renders 
the fibril rod-like, and the chiral nature of the beta strands 
restricts the number of filaments, 7, in a fibril [TtI]. In the 
fibril depicted below, 7 = 2. 
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(b) 
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figurational entropy for polymers. To simplify our the- 
oretical treatment, we will absorb these degrees of free- 
dom intrinsic to each monomer into and represent the 
monomer as a spherical particle. Since the free energy is 
defined up to addition of a constant, we will also set the 
monomer-solvent interaction energy to zero so that the 
free energy for free monomers becomes purely entropic. 

Without fibril-fibril interactions, the overall configura- 
tional partition function can be written as (cf. Appendix 
El): 



= n 



NJ. 



(1) 



where the prime in the product denotes the restriction 
that J2 S sN s = N, and 

Z s = exp[(s - 7)0 + s0] . 

To ease notation, we define two new parameters: 

X = a + + \n(u9) + ln(A" 3 *) (2) 
£ = 7 a + ln(u0) . (3) 

Namely, \ amounts to the sum of the monomer bind- 
ing energies through the A-type (first term) and P-type 
(second term) interactions, plus the entropic contribution 
(the third term) and the kinetic contribution (the fourth 
term); and £ amounts to the total longitudinal binding 
energies of the fibril, plus the entropic contribution. 

According to Eq. |T]), the Free Energy Density (FED), 
in the absence of fibril-fibril interactions, is expressed as: 



fo= dsn(s) [In n(s) - X s + C - 1] 



(4) 



where ksT is set to one and n(s) = N(s)/V with the 
unit volume set to be the volume of one monomer. Under 
this convention, n(s) is dimensionless and corresponds to 
the volume fraction. Note that in Eq. we have gone 
from a discrete description of the aggregation number to 
a continuous one. This assumption is valid if the mean 
aggregation number is large. 

To incorporate the steric interactions between fibrils, 
we employ the formalism developed in the study of cylin- 
drical micelles [H], H3|- Specifically, we model the free 
energy contribution of the fibril-fibril interactions as 



fir 



dsds'n(s)n(s')B(s,s') 



(5) 



where P(s, s') is the second virial coefficients of two rods 
of aggregation numbers s and s'. Denoting the diame- 
ter of the fibril by D and the length of a fibril with s 
monomers by L(s), we have 



B(s,s') = ^D 3 +^D 2 [L( S )+L( S ')]+DL{ S )L{s r )\sm< 



(6) 



where <f> is the angle between the two rods. Since the 
mean fibrillar length is much greater than D in our sys- 
tems of interest, we will ignore the first two terms in 
the second virial coefficients. Also, as a fibril is a linear 
structure, 



L{s) = (s 



(7) 



for some constant £. We shall from now on express the 
FED in terms of £• 

In the isotropic phase, the different directions of the 
rods are averaged over and so // is 



// = fo + 



irD( 2 



dsds' ss'n(s)n(s') , 



(8) 



In the nematic phase, the flexibility of the fibrils has to 
be taken into account to avoid length explosion in the 
nematic phase This can be done by incorporating 

the persistence length of the fibril, denoted by P, into 
the model. The resulting FED for the nematic phase is 
El: 



In 



P (s 
4A + 4A 



/jv = fo + J ds n(s) 

+ DC^/d 8 d S W„ W „( 8 <) , (9) 

where A is the deflection length of the fibril [28| . In par- 
ticular, A is related to the orientational order parameter, 
r], in the following manner [20l |: 



n ~ (l - 3A/P) 



(10) 



The length distribution that minimizes the above 
FEDs can now be found by using the Lagrange multi- 
plier method (e.g. see [H], EH)- For the isotropic phase, 
the distribution is : 



ni(s) = exp[-s/Si - £] 
Si = \fcif£ , 



(11) 
(12) 



where ci is the protein volume fraction and Si corre- 
sponds to the average aggregation number. Note that Si 
is above one only if cj > e - ^ , this therefore indicates that 
the CFC of the system is at e - ^ fl9j |. 

For the nematic phase, the distribution is 



4A* 

n N(s) = — cxp[-s/S N - £j 



S 



N 



Pc N ei 



4A* 



(13) 
(14) 



where, similar to the isotropic case, cn is the protein vol- 
ume fraction and Sm corresponds to the average aggrega- 
tion number. In the above equations, A* is determined by 
minimizing Eq. ([9]) with respect to A, and thus satisfies 
the following equation: 



(-N 



PA*e« 4A* 



Ccn , DC 2 



X*P 



■N 



= 



(15) 
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Due to the large magnitude of e^ in the systems that we 
are interested in (cf. Tabled, we will ignore the first term 
above and approximate A* as: 



A *3/2 = 



1 



2D(c N 



(16) 



Substituting Eqs CP and CE3]) into Eqs © and ©, the 
minimal FEDs for the two phases are: 



In 



■{xs 

4A* 



2Sj e" 4 + 



6 , ^DCcj 



4A* 



X ~ ^ ) Si + 2S N 



3 -£ 



-jv 



(17) 



(18) 



The FEDs above apply only to pure isotropic and ne- 
matic phases. To investigate the co-existence of the two 
phases, we denote by vi (vn) the proportion of isotropic 
(nematic) component in the system. The total FED is 
therefore 



/tot(c/,c/v) = «///(c/) +v N f N (c N ) , 

with the following conditions: 

UjCj + V N C N = c tot 
vi + v N = 1 . 

As vi, vn can be expressed in terms of cj, c/y: 



VI 



CN ~ Ctot 

cn - CI 



VN 



Ctot - C/ 
c N - Cf 



(19) 



(20) 
(21) 



(22) 



the total FED of the system is dependent only on cj 
and cat. The proportion of the isotropic/nematic com- 
ponent can now be obtained by minimizing the total FED 
with respect to ci and cn- This minimization problem 
does not admit an analytical expression and the graphi- 
cal method described in [2lJ is required (cf . Fig. [3]) . It is 
therefore worthwhile putting in experimentally relevant 
values into the model to reduce the number of variables 
to be analyzed. This leads us to the next section where 
the HL fibril solution is discussed. But before we do so, 
we note that as far as the length distribution and phase 
behavior are concerned, the energy strength, x, is irrele- 
vant. This is because the terms involving x i n the total 
FED, i.e., the first terms in Eqs (fT7)) and (fT8"|) , amounts 
to 



4A* q2 



(23) 



which is independent of c/ or cat. One consequence of 
this realization is that one will not be able to obtain a full 
picture of the overall fibrilization energy, which includes 
the longitudinal binding term, a, as well as the lateral 
binding term, /?, by studying the length distribution and 
the phase behavior of the system alone. 



FIG. 2: (Color online) Four beta sheet strands are shown with 
the fibrillar axis along the vertical direction. The hydrogen 
bonds are schematically displayed by the rods connecting the 
two pairs of beta sheets. On average, the amino-acids within 
a beta strand are about 0.35 nm apart. According to this 
structural model, each amino acid occupies on average (1 x 
0.35 x 0.48) = 0.17 nm 3 of space within the fibril. 



0.48 nm 




TABLE I: The various parameters involved in our model. The 
daggers indicate values predicted in this work. 



Properties 


Symbols 


Hen Lysozyme 


No. of AA 




129 


Persistence length 


P 


10 nm \8] 


Diameter 


D 


7.4 nm [10] 


Effective diameter 


-Doff 


13 nm 


Fibril length per monomer 


c 


0.5 nm 


Binding free energy 




20.6 f 


Lower cone, for phase sep. 


CA 


~ 0.65 mM [10] 


Upper cone, for phase sep. 


CB 


1.05 mM f 



III. HEN LYSOZYME AMYLOID FIBRILS 

We will now apply our theoretical formalism to a spe- 
cific system - the HL fibril solution. HL is a protein con- 
sisting of 129 Amino Acids (AA) and amyloid fibrils are 
observed to form when incubated at low pH and elevated 
temperatures [ToL |22| . The width and persistence length 
of a HL fibril is found to be 7.4 nm [l(| [22] and 10 iim [8[ , 
respectively. To complete the list of parameters involved 
in the model, we need to estimate £ (cf. Eq. ([7])) and £ (cf. 
Eq. |(3J|). Since the structural details of the HL fibrils are 
still lacking, we will estimate £ by employing the approx- 
imation adopted in in the study of elastic properties of 
amyloid fibrils. Specifically, each amino acid is assumed 
to occupy a volume of (1 x 0.48 x 0.35) = 0.17 nm 3 within 
the fibril. This assumption is motivated by the fact that a 
fibril constitutes mainly of cross-beta-sheet structure (cf. 
Fig. [2]). Given that D = 7.4 nm 10], the average length 
contribution to the fibril per monomer can be estimated 
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as follows: 



c 



129 x 0.17 nm 3 

(D/2)H 



0.5 nm. 



(24) 



In other words, each HL protein in the fibrils contributes 
on average 0.5 nm to the fibril's length. 

When the solution is at pH 2 in the presence of 100 mM 
NaCl, a HL protein carries a net positive charge of 19 [Icj |. 
Hence, a HL fibril has an average line charge density, v 1 of 
(19/C) = 38/nm. The electrostatic repulsion due to the 
fibrils' charge density can be accounted for by defining 
an effective diameter for the fibril, which is of the form 



D, 



off 



D 



In A + 0.577 + In 2- 1/2 



kD 



where 



Q 

A 



2N A e 2 I 

^2 



47reoefcsT 
7rL 2 



1 nm 



0.68 nm 



2Qne 



hD 



114 



(25) 

(26) 
(27) 
(28) 



In Eq. (|28[) . L ~ 300, and is estimated by solving the 
Poisson-Boltzmann equation for a charged cylinder in an 
ionic solvent (cf. Appendix [B| . Note that in the above 
equations, eo, e, Na, e, /, and Q are the vacuum 

permittivity, the dielectric permittivity of the solvent 
(taken to be 82 here), the Avogadro number, the ele- 
mentary charge, the ionic strength of the solvent in unit 
of mole/m 3 , the Debye screening length, and the Bjerrum 
length, respectively [HI, Substituting the values in 
Eqs fl2EJ| and (QSJ) to Eq. D cS can be calculated to 

be about 1.7 x D ~ 13 nm. We will employ this effective 
diameter in our FEDs shown in Eqs (fl"T|) and (|18p. 

To determine the only remaining parameter £, we will 
make use of the knowledge that the lower concentration 
for phase separation, ca, is measured to be around 0.6- 
0.7 mM Here, we set c A to be 0.65 mM for definite- 
ness. We then vary £ until c A obtained from the tangent 
method (illustrated in Fig. [3]) matches the assigned value 
of 0.65 mM. In doing so, we find that £ ~ 20.6 and the 
upper concentration for phase separation, cb, is about 
1.05 mM. We have now completely specified all of the 
model parameters, and based on these parameters, the 
variations of the various properties of the HL fibril so- 
lution with respect to protein concentration is shown in 
Fig.H 



IV. DISCUSSION 

Starting from a toy thermodynamic model , we have 
studied the isotropic-nematic phase transition in amyloid 
fibril solution by combining previous models for i) the 



FIG. 3: (Color online) The plot of /(c) + *c vs. c for the 
isotropic and nematic phases. The dotted line is tangent to 
both curves and the contact points (indicated by the arrows) 
correspond to the lower concentration, ca = 0.65 mM, and 
upper concentration, cb = 1.05 mM, for phase separation 



1.5 



x 10" 



f,(c)+Z c 
f N (c)+3C C 



0.5- 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1, 

c (mM) 



FIG. 4: The plots of a) nematic volume fraction, b) orien- 
tational order parameter, and c) average fibril length vs. HL 
concentration. 




nematic ordering in self-assembled linear structures, ii) 
the nematic ordering in charged rods, and iii) the elastic 
properties of amyloid fibrils. We then focus on HL fib- 
ril solution and estimated all of the parameters involved 
from experimental values. From the resulting parameter- 
free model, we deduced two main predictions: i) the up- 
per concentration for phase separation, cs, is 1.05 mM; 
and ii) the average fibril length varies with protein con- 
centration in the way depicted in Fig. Hfc) . In particular, 
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the average fibril length is predicted to be 1.21 /xm at the 
protein concentration of 0.5 mM. The first prediction on 
the upper concentration seems to be an underestimate 
of the upper concentration observed experimentally in 
(To| . This may be an outcome of further aggregation of 
fibrils. Indeed, it has been shown that amyloid fibril so- 
lutions tend to form gel at high concentration [12, HH . 
Gelation in fibril solutions points to possible attractive 
interactions between fibrils [24| , or contacts induced vis- 
coelasticity (25|. These effects are not captured in our 
FEDs, and may be the source of the discrepancy. On 
the prediction concerning the average fibril length, to the 
best of our knowledge, the average fibril length for HL 
fibrils has not been determined accurately and so this 
prediction remained to be verified. 

Besides protein concentration, salt concentration is 
also demonstrated to have a major effect on the onset of 
nematic ordering in [10( . Within our model, a decrease in 
ionic strength would increase the Debye screening length 
(cf. Eq. ([H])), an d hence also the effective diameter. This 
implies that a decrease in ionic strength would decrease 
the onset concentration for nematic ordering, which is 
indeed observed experimentally in flQ ] . 
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The first term corresponds to the kinetic contribution, 
the second term to the translational entropic contribu- 
tion with u being the roaming volume of each monomer 
within the fibril, the third term to the rotational entropic 
contribution with 9 being the roaming area on a unit 
sphere for the director of each monomer, and the fourth 
term to the binding energy, which is of the form: 

E s = (s- 7)0 + s(3 . (A3) 

In the above equation, a (/?) is the binding energy cor- 
responding to the A-type (B-type) interactions (cf. Fig. 
H}. Note also that the term s! in the denominator in Eq. 
(|A1|) disappeared in Eq. (|A2|) due to the fact that there 
are s! different ways of shuffling the monomers within the 
fibril. 

For a system with fibrils of variable lengths, we need 
to sum over all of the partition functions for the s-fibril, 
the total partition function is therefore: 

s 

where N s is the number of s-fibrils and the prime in the 
product denotes the restriction that J2 S sN s = N. 



The author thanks the Glasstone Trust (Oxford) and 
Jesus College (Oxford) for support. 



APPENDIX A: THE PARTITION FUNCTION 

The partition function for a fibril consisting of s 
monomers is 



Z„ = 



A 3s s! 



da:i • • • dx s dwi • • ■ dw s e~ 



-U{{x,w})/k B T 



(Al) 

where x (w) are the coordinates (directors) for the 
monomers, U({x,w}) is the potential function, and T s 
constrains the positions and directions of the monomers 
so that the aggregation is in the fibrillar form. Note 
that the prefactor corresponds to the kinetic part of 
the partition function such that A = hj \f%vmk^T 
is the dc Broglie thermal wavelength, and W = 
v/(27r) 5 (ksTYhhhlh? with 7, being the three principle 
moments of inertia [26( . 

In the mean-field limit where all monomer contribu- 
tions to the partition function are assumed to be identi- 
cal, the above integral can be partitioned into four terms 



x Vu 



8-1 



x AttO 



s-l x e E s 



(A2) 



translation rotation binding 



APPENDIX B: ESTIMATION FOR V 

The electrostatic potential ip for a long, cylindrical and 
charged rod in a solution with excess salt of ionic strength 
I is described by [27}: 



Id ( dtp\ _ 8irN A Ie 2 etp 
r dr \ dr J e$e ksT 



with the boundary conditions: 

ch/> 
dr 



4u 

r—D e eD 

lim ip(r) = , 



(Bl) 

(B2) 
(B3) 



where r denotes the distance away from the center line 
of the fibril. Note that we have assumed in Eq. (|F32|) that 
all of the charges are spread at the outer boundary of the 
cylinder. 

At a distance r far away from the rod (rer > nD/2 + 1), 



eij>(r) 
k R T 



TK (Kr) 



(B4) 



where Kq(.) is the modified Bessel function of the second 
kind The prefactor T can therefore be obtained by 
solving the differential equation Eq. (|B1[) 27]. 
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